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Abstract. Acoustic power and oscillation amplitudes of radial oscillations computed for a solar model are compared 
with solar seismic observations. The oscillations are assumed to be stochastically excited by turbulence. The 
numerical computation s are based upon a theo retical formulation of the power going into solar-like oscillation 
modes, as proposed by |5amadi fc Goupil (2001 ) in a companion paper. This formulation allows us to investigate 
several assumptions concerning properties of the stellar turbulence. We find that the entropy sour ce plays a 
dominant ro le in the stochastic excitation compared with the Reynold stress source, in agreement with 



Goldreich 



?t al. (1994). We consider several turbulent kinetic energy spectra suggested by different observations of the solar 



granulation. Differences between turbulent spectra manifest in large differences in the computed oscillation powers 
at high oscillation frequency. Two free parameters which are introduced in the description of the turbulence enter 
the expression for the acoustic power. These parameters are adjusted in order to fit to the solar observations 
of the surface velocity oscillations. The best fit is obtaine d with the kinetic energy spectrum deduced from the 
observations of the solar granulation by Nesis et al. (1993); the corresponding adjusted parameters are found to 



be compatible with the theoretical upper limit which can be set on these parameters. The adopted theoretical 
approach improves the agreement between solar seismic observations and numerical results. 
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1. Introduction 

It is currently believed that oscillations of solar-like stars 
can be excited by the turbulent convection of their outer 
layers. The physical description of this stochastic excita- 
tion which yields the power transferred into the oscilla- 
tions by the turbulence is not well known : the turbulent 



Reyn olds stress was first identified (Goldreich & Keeley 



1977] , hereafter GK) as the main source term. In the last 
decade, turbulent entropy fluctuations had been proposed 
as a possible additional source term ( Stein fc Nordlund 



1991; Balmforth, 1992 



Goldreich et al., 1994) 



In Paper I ( Samadi fc Goupil, 2001 ) we reconsider the 
excitation of stellar p-modes by turbulent convection. A 
new fo rmulation following the GK approach and taking 



Balmforth (1992| , hereafter B92) considered the corre- 
sponding contribution as negligible. However it was shown 
in Paper I that the entropy source term considered by 
GMK and B92 vanishes and a non-linear term involving 
the turbulent velocity field and the turbulent entropy rep- 
resents the actual entropy contribution. Inherent to an 
empirical description of the stellar turbulence, our theo- 
retical model for the acoustic power involves two free pa- 
rameters, the parameter A which is related to the rather 
arbritary definition of eddy correlation time and the pa- 
rameter f3 which relates the mixing length to the largest 
wavenumber in the inertial range. 

Several space seismic experiments are currently 
planed : COROT flBaglin fc The Corot Team 



ations contribution is proposed. In contrast with previ- 
ous works, the present formulation is valid for any turbu- 
lent time spectrum and turbulent energy spectrum. The 
entropy source term was considered as the main source 
term by Goldreich et al. (1994, hereafter GMK) while 



into a icount both Reynolds stress and entropy fluctu- jggjL MONS (|Kjeldsen fc Bedding, 1998Q and MOST 
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(Matthews, 1998). These projects will provide very accu- 
rate seismic data such as oscillation amplitudes and damp- 
ing rates for a large set of stellar targets. The oscillation 
powers deduced from these observations will provide valu- 
able constraints on the current physical understanding of 
the excitation by turbulent convection. In the meantime, 
we need to validate and constrain our theoretical approach 
with current solar observations. This is the primary goal 
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of this paper. We first compute the oscillation power for 
a solar model (Sect. 2) using conventional assumptions 
for the turbulent components and compare it with ob- 
servations. In Sect. 3 we investigate the influence of the 
atmosphere and of different turbulent components such as 
the time spectrum, the turbulent spectrum and the free 
parameters. 

In Sect. 4, we discuss on the best treatment for the 
turbulent components and then calibrate the correspond- 
ing formulation on solar observations i.e. determine the 
free parameters in view of applications to other solar-like 
stars. 

2. Oscillation power 

2.1. Theoretical formulation and observations 

The power going into each mode of frequency luq is ex- 
pressed as (see Paper I) 



p — r e 



(1) 



where T Wo = 2r\ is the linewidth of the mode, r\ is the 
damping rate and E Ma is the averaged mode energy given 

by 



Ew - ^ ( 
where 



M 



A \ 2 ) lUn 



1= / dmC4 
Jo 



(2) 



(3) 



is the mode inertia and (| A | 2 } is the mean-square am- 
plitude resulting from a balance between the excitation 
by the turbulent convection and damping processes. This 
amplitude can be expressed in a schematical form as 



oc rj 



M 



dm po w 



dr 



S(w, m) 



(4) 



where £ r is the radial displacement eigenfunction, po the 
density, w the vertical rms velocity of the convective ele- 
ments and § the turbulent source function including both 
the Reynolds and the entropy fluctuation contributions. 
Detailed expressions for (-*1 2 ) and 8 are given in Paper I. 
The source function involves the kinetic energy spectrum 
E(k), the entropy energy spectrum E s (k) and the time 
spectrum Xfc(w). The turbulent spectra in S are integrated 
over all eddy wavenumbers k, and 8 is in turn integrated 
in Eq. (|4j) over the stellar mass M. 

The time spectrum Xk{^) measures the temporal co- 
herence of eddies with given k. The linewidth of this func- 
tion is related to the time correlation Tk associated with 
eddies with wavenumber k. Tk is related to k and the eddy 
velocity Uk by (see Eq. 97 of Paper I) 



(5) 



where the free parameter A is introduced because the def- 
inition of the eddy correlation time and the evaluation 



of the eddy velocity are somewhat arbitrary. The eddy ve- 
locity Uk is given by Eq. (98) of Paper I. Several different 
time spectra are investigated later in this paper. In ad- 
dition, we evaluate effect of changes in the parameter A, 
which is usually set to one, for instance in B92. 

The kinetic energy spectrum E(k) describes the en- 
ergy carried by the turbulent eddies. In a highly turbu- 
lent medium, such as in the outer stellar convection zone, 
it is expected from turbulence theory that such a spec- 
trum exhibits an inertial range where the energy carried 
by the eddies with wavenumber k decreases as /c~ 5 / 3 in 
the Oboukhov-Kolmogorov theory. The wavenumber ko 
at which the turbulent cascade begins is not well deter- 
mined. This quantity has been related in Paper I to the 
mixing length A as 



(6) 



where p is a free parameter which has been introduced in 
order to gauge the arbitrary nature of such a relation. The 
injection region lies in the k < k range. Here we inves- 
tigate several turbulent kinetic energy spectra exhibiting 
different forms in the injection region. 

The mean-square amplitude is inversely proportional 
to the damping rate (e.g. Balmforth, 1992| ) and conse- 
quently the acoustic noise generation rate (or acoustic 
power) P Uo is independent of rj. Thus comparing the es- 
timated acoustic noise generation rates with observations 
avoids the additional uncertainties in the modeling of the 
damping rates (GMK) but requires accurate linewidth 
measurements and model computations. 

The mean-square surface velocity of a mode is given 
by the relation (see Paper I) 



1 



vi(uo) = ufe J (r a ) -(I A 



(7) 



where r s is the radius of the layer seen by the instru- 



ment. Observations of v a performed by Libbrecht (1988) 
are plotted in Fig. |l| with accompanying error bars. The 
surface velocity amplitude v s steeply increases with fre- 
quency at low frequencies, presents a maximum value 
around v ~ 3.2 mHz and undergoes a steep decline at high 
frequency. In Sect. 4 we make use of these three character- 
istics to constrain the turbulent ingredients and the free 
parameters entering the theory of the excitation process. 

2.2. Oscillation power in the solar case 

The solar model we consider has been calculated with 



the CESAM code (Morel, 1997) and appropriate input 
physics, described in details in Lebreton et al. (1999| ). In 



particular, convection is described according to the clas- 



sical mixing-length theory (Bohm - Vitense, 1958, here- 



after MLT) with a mixing- length A = a c H p , where H p is 
the pressure scale height and a c is the mixing-length pa- 
rameter. The external boundary conditions are obtained 
through a T(t) law d erived from ATLAS9 model atmo- 
spheres ( Kurucz, 1991 ). The solar metal to hydrogen mass 
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ratio (Z/X)q = 0.0245 comes from Grevesse & Noels 
(1993j)] The calibration of the solar model, in luminosity 
(L Q = 3.846 1 33 erg/s) and radius (R Q = 6.9599 10 10 cm) 
for an age of 4.6 Gyr, fixes the initial helium content 
Y = 0.2682, metallicity Z = 0.0175 and the MLT pa- 
rameter a c = 1.785. 

The estimated cut-off frequency is about v c ~ 
5.32 mHz. The oscillation eigenfunctions are obtained with 



Leon 



the adiabatic pulsation FILOU code of Tran Minh , 
(19950 

We evaluate the root mean-square surface velocity v s 
in Eq.(Q) where the oscillation power P Wo is computed 
accordi ng to Eqs(|l] , p| , |) and rj is given by solar observa- 
tions of |Libbrccht (1988 ). These solar observations sample 
the solar atmosphere at the optical depth T5000 — 0.05. In 
our solar model this optical depth is localized approxima- 
tively at 140 km above the solar photosphere and v s is 
thus computed at this depth. 

Fig. [l shows the resulting values of v s assuming the 
Kolmogorov spectrum as defined in Eq. (j|) . 

As in GMK and B92, it is assumed that the entropy 
spectrum exhibits the same behavior than the kinetic en- 
ergy spectrum. Thus the entropy energy spectrum E s (k) 
considered in the computation is at this stage simply pro- 
portional to the kinetic energy spectrum E(k) (Eq. 107 
of Paper I). The parameter A is first set to one and the 
parameter j3 is set to 1.9, according to the theoretical es- 
timate obtained in Paper I. The anisotropic factor $ in- 
volved in the formulation may be considered as a free pa- 
rameter. However in order to be consistent with Bohm- 
Vitense's MLT considered here, we set $ = 2 in the 
present work. $ affects the v dependence of the power 
in a similar way as A and /3. Changing its value is there- 
fore equivalent to change values of A and (3, which will be 
adjusted later (section ||). 

At low frequencies, v s increases because the inertia 
of the modes decreases with frequency and because the 
derivatives of the displacement eigenfunctions increase 
with frequency (see GMK). At high frequency, v s de- 
creases because the depth of the excitation region and the 
efficiency of the excitation decrease with frequency. This 
can be explained as follows: let be the characteristic life- 
time of an eddy of wavenumber fc. The major contribution 
to mode excitation comes from the eddies with Tktoo < 1 
( |Goldrcich fc Kumar, 1990| , Paper I). When ujq increases, 
the region where r^wo ^ 1 becomes thinner and is confined 
at the top of the convection zone. Futhermore, contribu- 
tive eddies have smaller size and are less energetic. Thus 
P Ua becomes smaller as loq increases. 

Including only the Reynolds stress leads to amplitudes 
much smaller than those observed. The contribution of 
the entropy source term is found to be larger than the 
Reynolds stress contribution and increases the amplitudes. 
This result is in agreement w ith the conclusion of GMK 
and with the simulations of Stein fc Nordlund (1991 ). 




Fig. 1. Root-mean-square velocity v s (u>o): dots with ac- 
companying error bars represent the solar measurements 
of Libbrecht (1988) ) . The continuous curve shows velocity 
amplitudes computed according to Eq.(0) using both the 
Reynolds and the entropy contributions and the dashed 
curve using only the Reynolds stress contribution. A 
Kolmogorov energy spectrum is assumed. Damping rates 
are deduced from Libbrecht (1988) and we assumed A = 1 
and 13 = 1.9. 



the observed ones and large differences are observed at 
high frequency in terms of frequency dependence between 
the amplitudes computed using the Kolmogorov spectrum 
and the observations. Also, the maximum amplitude peaks 
at a frequency which is smaller by about 0.1 mHz, com- 
pared to the observed position of the maximum. 

3. Impact of different assumptions 

The above calculations were performed considering stan- 
dard assumptions for the turbulent spectra and the free 
parameter A. We now investigate effects of changing the 
turbulent properties, the free parameters A and (3, the 
traitement of the atmosphere and of the depth in the 
atmosphere at which the expected surface velocity v s is 
computed. 

3.1. Turbulent kinetic energy spectrum 

In Sect. 2 we considered the Kolmogorov Spectrum (KS 
hereafter) 



E(K) 



11, 



K- 5/3 for K > 1 
fc 



(8) 



However, some discrepancy still remains between the the- 
oretical estimates and the observations. In particular the 
theoretical absolute amplitudes remain much smaller than 



where K = fc/fco, fco = 2tt/(/3A) and Uq is proportional to 
the mean square turbulent velocity (see Paper I). 

Observations of solar granulation show a change of 
slope of the kinetic energy spectrum between the inertial 
range where the spectrum obeys the Kolmogorov law and 
the large scale range (small wavenumbers fc) correspond- 
ing to the energy injection region. 
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Fig. 2. The RKS, KS, SS, BKS and NKS versus the nor- 
malized wavenumber K. The NKS has been shifted up in 
order to distinguish it from the BKS for K > 1 

In order to take large scale eddies into account, we 
consider first the Spie gel spectrum (S S hereafter) for com- 
parison with B92 and Houdek (1996|). 



Fig. 3. Amplitudes computed with the Reynolds stress 
contribution only and for the kinetic energy spectra of 
Fig. |^. A gaussian time spectrum has been adopted. The 
curves are normalized to their maximum value and have 
intentionally been shifted by +0.1 mHz in order to match 
the position of the observed maximum. All other assump- 
tions are those of Fig. [l| 



E(K) =18^ (2Ky 5/3 J(2K) for K > 0.5 
k 

with 
J(K) 



(9) 



which mimics the solar kinetic spectrum observed by Nesi: 
et al. (1993Q . 



1 _ - -(1 - K- 1 ^) 



(10) 



This spectrum takes into account the non-linear in- 
teraction between turbulent modes of low wavenumber 
( [Bpiegcl, 1962 ; Stein, 1967) and is built for very large 
Rayleigh numbers. 

In order to match the high- resolution observa- 
tions o f the solar granulati on by Roudier fc Mullcr 
(19860 |Musielak et al. (1994j ) has proposed the "Raised 
Kolmogorov Spectrum" (RKS hereafter) 



E(K) 



0.54 -0^- 5 / 3 for K > 1 
k 



(11) 



0.54 -°K-°- 7 for 0.2 < K < 1 



and as a possible alternative the "Broad Kolmogorov 
Spectrum" (BKS hereafter) 



0.652 J2. K~ 5 ' 3 for K > 1 



E(K) 



fe 



0.652^° 

k 



(12) 



for 0.2 < K < 1 



We also consider the "Nesis Kolmogorov Spectrum" 
(NKS, hereafter) 



E(K) 



0.655 ^" 5/3 for K > 1 
u 2 

0.655 for 0.7 < K < 1 



(13) 



The above five kinetic energy spectra are plotted in 
Fig. ||. All these spectra exhibit the same —5/3 slope at 
large wavenumbers. The RKS, the BKS, the SS and the 
NKS take into account large-scale eddies lying in the in- 
jection region. 



As stressed by Rieutord et al. (2000 ) , differences at low 
wavenumber between the kinetic spectra obtained from 
observations of the solar granulation are a consequence of 
uncontrolled data-averaging procedures. Properties of the 
kinetic spectrum at low wavenumber (mesogranulation) 
are thus not yet well depicted by the observations of the 
solar surface. In the present paper we nevertheless study 
the different spectra of Eqs.(|ll-13) in order to gauge the 
sensitivity of the acoustic power emission to the properties 
of the turbulent spectrum in the injection region. 

Fig. H presents amplitudes computed with the sole 
Reynolds stress contribution and for the different kinetic 
energy spectra of Fig. ||. Only small differences in the 
velocity amplitudes are observed arising from different 
spectra at low frequencies. The main differences are ob- 
served at high frequency where the RKS and BKS exhibit 
a steeper slope than the KS and SS. These features are 
explained below. 

The velocity amplitudes obtained by B92 using the KS 
are shifted by a larger amount toward lower frequencies 
compared to ours and diff erences betwee n the KS and SS 
are found to be larger in Houdek (1996 ), in which com- 
putations are based on the formulation of B92. The last 
feature is due to the way B92 implemented the SS and the 
first feature is due to the way we normalized the turbulent 
energy spectra (Eq. 63 in Paper I). 
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Fig. 4. Velocity amplitudes computed with both the 
Reynolds and the entropy contributions and for different 
kinetic energy spectra. A gaussian time spectrum has been 
adopted. The curves are normalized to their maximum 
value and have intentionally been shifted by +0.1 mHz 
in order to match the position of the observed maximum. 
Dots are observational points as defined in Fig.|l|. The free 
parameters are set to A = 1, (3 = 1.9. 

With our formulation, velocity amplitudes using the 
SS are similar to those using the KS. Therefore we no 
longer consider the SS in the remaining part of this work. 

From here on, amplitudes calculations include both 
Reynolds stress and entropy contributions. Fig. ^ shows 
amplitudes computed with the free parameters set to 
A = 1, j3 = 1.9 with <I> = 2. Large differences are found 
at high frequencies arising from the use of different tur- 
bulent spectra. The RKS leads to larger amplitudes and 
fits better the solar observations at high frequencies. At 
low frequencies, the RKS overestimates v s more than the 
estimates including only the Reynolds stress contribu- 
tion. The RKS contains a significant amount of kinetic 
energy at low wavenumbers. At low frequencies the ex- 
citation process is mainly due to the large scale eddies. 
Furthermore it is stressed in Paper 1 that the entropy 
contribution is more efficient for long-period oscillation 
modes. Therefore the RKS injects through the entropy 
contribution a larger amount of energy to the low fre- 
quency modes than the other spectra. For the other spec- 
tra we found few differences in terms of frequency depen- 
dence between computations including only the Reynolds 
stress contribution and including both Reynolds stress and 
entropy contributions. 

The computed P uo is plotted in Fig. ^ in order to em- 
phasize the changes due to the use of different turbulent 
spectra. Also shown is the power P UJo derived from the 
observations according to Eq.(^) where the radius r s cor- 
responds to the layer at which v s is measured. We observe 
large differences between the oscillation spectra P Uo . At 
high frequency, the KS strongly overestimates P Wo while 
the RKS and the BKS are closer to the oscillation power 
derived from the observations. Thus the RKS and BKS 




v (mHz) 



Fig. 5. P uo for different kinetic energy spectra are pre- 
sented in log — log representation. Dots represent the os- 
cillation power derived from the solar observations with 
the help of Eq. (R). Corresponding errors bars are plotted. 



inject relatively less energy into the high frequency oscil- 
lations. This is explained as follows: contributive eddies 
are those for which t^cjo J$ 1. At a given wavenumber k 
the RKS induces a smaller velocity Uk for the considered 
eddy. As w (fcufc) -1 for a given oscillation frequency, 
the RKS then involves eddies with larger wavenumbers 
than the KS. Eddies with a high wavenumber carry less 
energy and consequently the RKS injects relatively less 
power into the oscillations of high frequencies than the 
KS does. 

At low frequencies smaller differences are observed be- 
tween the spectra because for small values of the oscilla- 
tions frequencies all the eddies of the spectra are involved 
in the excitation process. However, the RKS overestimates 
P Wo while the BKS leads to an intermediate result between 
the RKS and the KS. 

Fig. H compares the sizes of the excitation region when 
using the KS and the RKS for several oscillation modes. It 
shows that the depth dependence of P(u>) behaves accord- 
ing to the depth dependence of the associated displace- 
ment eigenfunction. For higher radial order, the mode is 
more concentrated toward the surface and consequently 
the excitation region is shallower. The extension of the 
excitation region is about twice the size with the KS 
than with the RKS; i.e. « 1.5 Mm with the KS against 
ps 0.8 Mm with the RKS and for a mode of order n = 20 
and a frequency i/«3 mHz. This can be explained as fol- 
lows: For a given wavenumber k, the correlation time of 
an eddy, is larger with the RKS than with the KS because 
the RKS induces smaller velocity Uk at a given wavenum- 
ber k. Therefore the region where TkLOQ < 1 is smaller with 
the RKS. This explains the strong differences in terms of 
depth of excitation found between the KS and the RKS 
(Fig. ^). The BKS induces an intermediate depth of exci- 
tation between the RKS and the KS (not shown here). 
Therefore the RKS induces a thinner excitation region 
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Fig. 6. Depth dependence of the injected power in the solar outer layers computed for the oscillation modes of radial 
order n = [10, 15, 20, 25, 30]. The abscissae are the stellar radius from center in Mm. The oscillation power is plotted 
for the KS (top) and the RKS (bottom) with A = 1 and (3 — 1.9. All other assumptions are the same as in Fig. 4. 



which is c loser to the one found i n the simulation per- 
formed by Stein fc Nordlund (2000 ) . 



3.2. Turbulent entropy spectrum 



It has been assumed above that the entropy spectrum ex- 
hibits the same behavior as the kinetic energy spectrum. 
However, in paper I, we stress that the entropy spectrum 
exhibits the same behavior as the tempera ture spectrum 
and the turbulence theory ( Lcsieur, 1997 , Chap VI-10) 
predicts a specific spectrum for the turbulent tempera- 
ture fluctuations which derives from the scalar nature of 
the temperature. We therefore proposed the "conductive" 
entropy spectrum derived from the theoretical spectrum 



expected for the temperature fluctuations (Lesieur, 1997) 



a s 2 



E,(K) 



aps 2 
u^ko 



E{K) 

Kc 
K 



for 1 < K < K c 



E(K) for K c <K 



(14) 



where clq is a normalization factor, E{K) is the kinetic 
energy spectrum, s 2 is the entropy scalar variance (Eq. 88 
of Paper I) and Kc = k c /ko where k c is the wavenum- 
ber which separates the inertial-convective range from the 
inertial-conductive range. When the kinetic energy spec- 
trum E(k) obeys the Kolmogorov scaling law, the entropy 
spectrum verifies 



E S (K) 



aps 2 

kp 
aps 2 



#-5/3 



for I < K < K, 



c 



K% K~ 17/3 for K c < K 



(15) 
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Fig. 7. Amplitudes computed with the "Kolmogorov 
Conductive spectrum" (KCS) for different values of Kq- 
All other assumptions are those of Fig. 4. 



This spectrum, that we call the "Kolmogorov Conductive 
spectrum" (KCS, hereafter), matches observations of in- 
tensity fluctuations of the s olar granulation p erformed by 
Espagnet et al. (1993|) and |Nesis et al. (1993| ). 

There is no theoretical constraints on the value of Kc ; 
therefore Kc may be considered as a free parameter of our 
formulation. We thus investigate different values of Kc 
with the use of the KCS (Eq.15). The results are presented 
in Fig. [| 

We observe no change in the shape when varying Kc, 
only a decrease of the global amplitudes with increasing 
Kc- As Kc decreases, the major part of the energy in the 
entropy energy spectrum is concentrated at low wavenum- 
bers. Eddies with small wavenumbers carry more energy 
than those with larger wavenumbers. Therefore a decreas- 
ing value of Kc induces an increase of P Uo and thus the 
shape of v s (loo) remains unchanged. We found similar re- 
sults with the entropy "conductive" spectra based on the 
RKS, BKS and NKS (not shown here) and derived from 
Eq.(g). 

Therefore we cannot constrain the value of Kc from 



the seismic solar observations. However Espagnet et al 
(1993|J] shows that the incrtial-convective domain of the 
temperature spectrum starts from granular scales of about 
ko « 3M _1 m while the inertial-conductive range starts 
from k c w 6 M _1 m. This leads to a value Kc = 2. 

In the remaining part of this work, the computation of 
P Uo will be performed assuming the "conductive" spectra 
for the entropy contribution with K c = 2. 



3.3. Time spectrum 

The time spectrum function of eddies with given k is usu- 
ally modelled by a gaussian function in which lincwidth is 



proportional to the time correlatio n of the eddy ( poldrcich 
k Kecley, 1977j |Balmforth, 1992|, Paper I). As proposed 



sider below other forms for the time spectrum Xk(u) such 
as the modified gaussian form 

(16) 



\/7T Wi 



and an exponential form 
1 



Xfc(<*>) 



-\u/u k \ 



2wj 



(17) 



The quantity lj^ controls the width of the eddy time spec- 
trum and is therefore related to the eddy correlation time 
as Wfc = 2/rfe (see Paper I). All time spectra are normal- 
ized such that 

f + OO 

Xk {uj)duo = l. (18) 

This choice of normalization differs from the one consid- 
ered by Musielak et al. (1994), but is in aggrement with 



by [Btein (19671) and |Musiclak et al. (1994| ) we also con 



Tennekes fc Lumley (1982j ) requirements . 

Fig. |8| shows the normalized velocity amplitudes v s 
computed with the RKS and for different time spectra 
: the gaussian form, the modified gaussian and the ex- 
ponential forms. Both the exponential and the modified 
gaussian forms yield amplitudes which are too large com- 
pared with the observations at high frequencies while the 
gaussian form gives rise to amplitudes which are closer 
to the observed amplitudes. At low frequencies both the 
modified gaussian and the exponential forms fall slightly 
closer to the observations than the gaussian form. In con- 
trast, computed v s assuming the KS leads to the opposite 
behavior for the gaussian form, i.e. at low frequencies, 
computed v s with the gaussian form are closer to the ob- 
servations than the other forms (not shown here). 

The strongest disagreement between observed and 
computed surface velocity amplitudes using KS and RKS 
is obtained at high frequency for the exponential form. 
Therefore we discard this time spectrum from now on. We 
cannot clearly discriminate between the modified gaussian 
form and the gaussian form and choose to adopt the gaus- 
sian form. 

3.4. Influence of the atmosphere 

We have computed two solar models, one with the 
Eddington classical gray atmosphere and one with the 
Kurucz ATLAS 9 (Kurucz, 1991) solar model atmosphere. 
The differences in the computed powers (not shown here) 
are found to be much smaller than effects due to the choice 
of the kinetic spectrum and than the errors in the obser- 
vations. 

On the other hand, the measured velocity amplitudes 
v s depend on the height in the atmosphere where the ob- 
servations are sampled. Fig. |] shows theoretical velocity 
amplitudes calculated at different heights in the atmo- 
sphere. 

Changes in the eigenfunction profiles at the solar sur- 
face are larger with high frequency modes. Therefore in- 
creasing the height in the atmosphere where v s is com- 
puted leads to larger changes at high frequencies (y > 
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gaussian 
modified gaussian 
exponential 




Fig. 8. Computed v s (cuq) assuming the RKS is plotted for 
different frequency factors Xk(^)- All other assumptions 
are those of Fig. | 





Qx = 2.9 cm/s 

0K = 3.0 cm/s 

ox = 3.2 cm/s 

ay = 3.4 cm/s 

x = 3.7 cm/s 



[mHz] 



Fig. 9. Velocity amplitudes v s computed at different levels 
in the atmosphere and using the RKS. Other assumptions 
are those of Fig. 4. 



3.5 mHz). Almost no effect is observed at low frequencies. 
This shows that the physics of the excitation process is 
more easily constrained by the observations at low fre- 
quencies (y < 3.5 mHz). 



the parameter A. In terms of frequency dependence, this 
is equivalent to fixing A and to varying the parameter [3. 

Fig. [To] shows velocity amplitudes v s computed when 
assuming the KS, the RKS and the NKS for different val- 
ues of A. The value of A (or (3) controls the extension of 
the excitation region through Tjt . Small values of A lead to 
a larger excitation region and therefore favor long-period 
oscillation modes. This explains the larger frequency shift 
of the maximum amplitude between the computations and 
the observations obtained with A < 1. Changing A has 
smaller effects at low frequencies than at high frequencies 
because at low frequencies all the eddies in the spectra 
are involved in the excitation process and the frequency 
dependence is mainly controlled by the mode inertia and 
the derivatives of the eigenfunctions. Changes are never- 
theless larger when using the RKS than the NKS and the 
KS. 

For A = 0.5 and using the KS, the frequency shift be- 
tween the numerical results and the observations is about 
0.25 mHz and the computations overestimate the ampli- 
tudes at low frequencies. In that case, the velocity v s {ojq) 
is close that found by B92 and Houdek (1996) which as- 
sumed A = 1. This result is consistent with the Osak 



(1993) statement that the B92 formulation leads to a 
source of excitation that is too extended. Our estimates 
with A = I and using the KS are closer to the observa - 
tions than those performed by B92 and Houdek (1996 ). 
This difference is due to the way we have normalized the 
turbulent energy spectra. Therefore our formulation leads 
to a smaller depth of excitation and is in better agreement 
with the observations. 

We observe large changes at high frequencies with A 
for all the spectra. Increasing A induces a steeper slope 
at high frequency. This can be explained as follows: con- 
tributive eddies are those for which r^wo < I. These eddies 
are the eddies for which k > fco(A/3 ujq/ KuKtotCf^ 2 ■ As A 
increases, only eddies with higher wavenumbers k are in- 
volved in the excitation process. Therefore, because eddies 
with high wavenumbers carry less energy, P UJo is smaller 
with increasing value of A. 

The oscillation amplitude increases as /3 4 , hence (3 
plays a major role in the control of the maximum am- 
plitude. A theoretical constraint can be obtained for (3 as 
(3X < 2.7 $ 1/2 (see paper I). For $ = 2, we have j3\ < 4. 



3.5. Effects of the free parameters A and (3 

With fco given by Eq.6, the eddy time correlation length 
(Eq.5) can be related to j3 and A as 



Tk 



\(3 



uj\Ku 



K 



(19) 



where wa = 2ttuq/ A and uk — Uk/uo. Thus both (3 and A 
are involved in the expression of the amplitude (Eq. 101, 
Eq. 102 and Eq. 103 of Paper I). The product (3 A controls 
the frequency-dependence of the oscillation power. Fixing 
the value of j3 to 1, we investigate the effects of changing 



4. Solar calibration 

We need to determine the free parameters A and (3 by 
calibrating the theoretical amplitudes on the solar obser- 
vational ones. For each of the various kinetic and "con- 
ductive" spectra, we first adjust (3X in order to match as 
closely as possible the position of the maximum and when 
possible the frequency dependence of vJujq) to the ob- 



served surface velocity given by Libbrccht (1988 ). We next 
adjust (3 and A in such a way as to obtain the observed 
maximum value of 18cms _1 . Table |l| gives the values of 
the adjusted parameters. 
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. [mHz] 



Fig. 10. Surface velocity v s (u>o) computed with j3 = 1 for 
different values of the parameter A. The top panel concerns 
the KS, the middle one the RKS and the bottom one the 
NKS. All other assumptions are those of Fig. 0. 



Fig. |ll| shows amplitudes computed with the RKS, KS, 
BKS and NKS with (3 and A set to their adjusted values. 
The KS and NKS yield much better fits to the seismic 
observations than the other spectra. 



Observations of solar granulation by Espagnet et al 
(19930show that fco — 3 Mm 1 whereas a value of fcp ~ 
7 Mm -1 can be inferred from the observations of Ncsis 



et al. 

(Oil)- 

These values correspond to (3 ~ 8 and (3 ~ 4 
respectively. The values of (3 adjusted with the RKS, the 
BKS and the NKS are closer to the value of (3 suggested 



Table 1. Values of the adjusted parameters (3 and A. The 
wavenumbers fco of the largest eddies in the incrtial range 
have been evaluated at the solar radius according to the 
relation kg = 2ir/[3A. We have assumed the "conductive" 
spectra and Kc — 2. 



spectrum 


px 




A 


fco (Mm" 1 ) 


RKS 


1 


3.53 


0.28 


7.3 


BKS 


1 


4.44 


0.22 


5.8 


NKS 


3.8 


5.00 


0.76 


5.1 


KS 


5.7 


9.51 


0.60 


2.7 




Fig. 11. Velocity mplitudes computed with the RKS, 
BKS, NKS and KS with values of the parameters (3 and 
A as given in Table 1. [j]. All other assumptions are those 
of Fig. § 



by the observations of Nesis et al. (1993) while the value 



suggested by the observations of Espagnet et al. (1993 ) 
favors the KS. Therefore (3 is not well constrained by the 
observations of the solar granulation and this cannot help 
to discriminate between the different turbulent spectra. 

The best fits are obtained with the NKS and the KS. 
The power computed with the KS seems to fit best the 
observations; this can be explained by the existence of 
compensating effects. We indeed know from the observa- 
tions of the solar granulation that kinetic energy is carried 
by large-scale eddies. This corresponds to the k— injection 
region lying at scales above the inertial range, which is not 
described by KS. Furthermore, the theoretical constraint 



imposed on (3\ (see Sect. 3.5) suggests that the adjust- 
ment value of (3\ found by fitting the amplitudes with the 
KS to the observations is not valid, while the other spec- 
tra verify this constraint. One is therefore left with the 
NKS, A = 0.76 and (3 = 5.0 as our best fit of the solar 
data. Howev er, the fit is not per fect at low frequencies. 
We show in Samadi et al. (2001 ) that this is related to 
the treatment of the static model and the corresponding 
oscillations. 
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5. Conclusion 

Acoustic power and oscillation amplitudes of radial os- 
cillations computed for a solar model are compared with 
solar seismic observations. With the adopted approach for 
computing the acoustic power (paper I), we find that the 
entropy source term significantly contributes to the exci- 
tation process of solar-like oscillations by turbulent con- 
vection, in agreement with the results of GMK and simu- 
lations by |Stein fc Nordhmd (199l|). 



Our generalized formulation for the oscillations power 
Pujo injected into the modes allows us to investigate effects 
of changes in the turbulent components. 

In the solar case, several kinetic energy spectra have 
been considered: the Kolmogorov spectrum (KS); the 
"Raised Kolmogorov Spectrum" (RKS) which is derived 
from observations of the solar granulation by Eloudier & 
Muller | (19861 ); the "Nesis Kolmogorov Spectrum" ( NKS) 
which is derived from observations performed by Nesis 
et al. ( ,993 ). We have, in addition, considered the "Spiegel 
Spectrum" (SS) which takes into account the non-linear 
interaction between turbulent modes of low wavenumber 
and the "Broad Kolmogorov Spectrum" (BKS) which ex- 
hibits intermediate properties between the SS and the 
NKS at low wavenumbers. 

In a first step, the value of the free parameter (3 is 
set to a value loosely determined using theoretical argu- 
ments and A = 1 is assumed as in B92. The oscillation 
power spectra computed with the KS, BKS, NKS and 



RKS significantly differ. In contrast with Houdek's (1996) 
results, only small differences between the KS and SS are 
found. At high frequency, the RKS gives rise to a frequency 
dependence of the velocity which more closely fits the 
frequency dependence of the solar observations than the 
other spectra and in addition induces larger amplitudes. 
Moreover, the RKS leads to a depth of excitation of th e 
same order as the one found by |Stcin fc Nordhmd (2000| ). 
This depth of excitation is found to be much smaller than 
the one assoc iated with the KS. This result is agreement 
with Osaki's ( 1993| ) statement that the source of excitation 
should be more concentrated than the extended source dis- 
tribution found by B92. 

As a model for the time spectrum of the turbulent ed- 
dies, we find that the gaussian form, which is usually con- 
sidere d, and the modified gaussian proposed by Musielak 
et al. (1994) , both agree with the observations while the 
exponential form - also proposed by Musielak et al. (1994) 
- is not satisfactory. 

We have assumed an entropy spectrum which is ex- 
pected from the turbulent theory. This spectrum presents 
a steeper slope in the inertial-conductive region, i.e. above 
the wavenumber k c . Only changes in the global amplitudes 
are observed when changing k c : the amplitude increases 
with decreasing values of k c . Thus the value of k c is diffi- 
cult to constrain because the global amplitude depends on 
the treatment of convection. Therefore the entropy spec- 
trum expected from the theory cannot be validated with 
the solar seismic observations. The value of k r has thus 



been set to the value suggested by the observations of the 
solar granulation by Espagnet et al. (1992). 

Assuming a theoretical value for the free parameter (3 
(as suggested in Paper I) for each spectrum, we find that 
the absolute amplitude remains much smaller than the ob- 
servations (about 3—5 cm/s). Because the formulation is 
very sensitive to /?, the observed differences in the maxi- 
mum amplitude between the theoretical estimates and the 
observations are mainly related to the uncertainties in the 
determination of (3. Observations of solar granulation are 
not accurate enough to provide a very precise value of 
the parameter (3. Also the product /3A controls the depth 
of excitation and therefore also the frequency dependence 
of the expected power. Accordingly, (3\ is given a value 
for each spectrum which best matches the frequency de- 
pendence of the power spectrum at low frequencies to the 
observations. We next determine the value of (3 and A in 
order to adjust the maximum amplitude to the observa- 
tions. This in turn defines the eddy time correlation and 
the scale at which energy injection is initiated. 

Doing so, all the turbulent spectra give rise to a sim- 
ilar behavior for P^ at low frequencies but the KS and 
NKS fit much better the observations at high frequencies. 
The values of (3 and A resulting from the fitting performed 
with the RKS, BKS and NKS are found to be compatible 



with the Nesis et al. (1993) observations of the solar gran- 
ulation. For the KS we find much larger values of (3 and 
A. The corresponding value of (3\ is not compatible with 
the estimated upper limit value of (3X. This suggests that 
the NKS is in better agreement with the observed surface 
velocity spectrum of the Sun 

The agreement in using the NKS between theoretical 
results and solar observations is found to be much bet- 
ter than that obtained in previous work. In particular at 
high frequencies, the power derived from the observations 
exhibits a slope of about —6 above v = 4.5 mHz. 

GMK predicts a slope of about —4.4 and a slope of 
about —1.0 is found with B92's formulation. In the present 
work, we find a slope \og(P UJo ) vs log(^) of about —6 in 
the frequency domain v > 4.5 mHz. 

It is not obvious that the above conclusions hold true 
for other solar-like oscillating stars. In particular one ex- 
pects Puj to exhibit different properties for stars with ro- 
tation and magnetic activity which are different to that 
of the Sun and the NKS might well not not suitable to 
represent the actual turbulent spectrum in these stars. 

Assuming the same properties of turbulence i.e. same 
values of the free parameters for a given turbulent spec- 
trum, it is possible to apply our calibrated formulation 
to several potentially solar-like oscillating stars in order 
to investigate the consequences of changing the turbulent 
ingredients and to establish constrains on the stellar tur- 
bulent spectrum ( jBamadi et al., 2000| , |2001 ). 

On the theoretical side, further tests and validation are 
needed. In particular, as suggested in Paper I, 3D simula- 
tions of the stellar convection will enable us to test some of 
the assumptions and approximations entering the present 
formulation (Samadi, 2000). In addition, such simulations 
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may enable us to constrain more accurately the values of 
the free parameters than has been achieved here with solar 
seismic data. 

Acknowledgements. We thank F. Tran Minh and L. Leon for 
the use of the FILOU pulsation code and we are indebted to 
J. -P. Zahn for useful discussions. 



Stein, R. F. & Nordlund, A. 2000, ApJ(in press) 
Tennekes, H. & Lumley, J. 1982, A first course in turbu- 
lence, 8 th edn. (MIT Press) 
Tran Minh, F. & Leon, L. 1995, in Physical Process in 
Astrophysics, 219 



References 

Baglin, A. & The Corot Team. 1998, in IAU Symp. 185: 
New Eyes to See Inside the Sun and Stars, Vol. 185, 301 
Balmforth, N. J. 1992, MNRAS, 255, 639 
B6hm - Vitcnsc, E. 1958, Zeitschr. Astrophys., 46, 108 
Espagnet, O., Mullcr, R., Roudicr, T., & Mcin, N. 1993, 

A&A, 271, 589 
Goldrcich, P. & Keeley, D. A. 1977, ApJ, 212, 243 
Goldreich, P. & Kumar, P. 1990, ApJ, 363, 694 
Goldrcich, P., Murray, N, & Kumar, P. 1994, ApJ, 424, 
466 

Grevesse, N. & Noels, A. 1993, in Origin and Evolution 
of the Elements, ed. N. Prantzos, E. Vangioni-Flam, & 
M. Casse (Cambridge University Press) 

Houdck, G. 1996, PhD thesis, Institut fur Astronomic , 
Wien 

Kjeldsen, H. & Bedding, T., cds. 1998, The First MONS 

Workshop: Science with a Small Space Telescope, ed. 

H. Kjeldsen & T. Bedding, Aarhus Universitet. 
Kurucz, R. L. 1991, in Stellar Atmospheres - Beyond 

Classical Models, 441 
Lebreton, Y., Perrin, M., Cayrel, R., Baglin, A., & 

Fernandes, J. 1999, A&A, 350, 587 
Lcsieur, M. 1997, Turbulence in fluids (Kluwer Academic 

Publishers) 
Libbrecht, K. G. 1988, ApJ, 334, 510 
Matthews, J. M. 1998, in Structure and Dynamics of the 

Interior of the Sun and Sun-like Stars, 395 
Morel, P. 1997, A&AS, 124, 597 

Musiclak, Z. E., Rosner, R., Stein, R. F., & Ulmschneider, 

P. 1994, ApJ, 423, 474 
Nesis, A., Hanslmeier, A., Hammer, R., Komm, R., 

Mattig, W., & Staiger, J. 1993, A&A, 279, 599 
Osaki, Y. 1993, in ASP Conf. Ser. 40: IAU Colloq. 137: 

Inside the Stars, 512-520 
Rieutord, M., Roudier, T., Malherbe, J. M., & Rincon, F. 

2000, A&A, 357, 1063 
Roudier, T. & Mullcr, R. 1986, Solar Phys., 107, 11 
Samadi, R. 2000, PhD thesis, Universite Paris 6 
Samadi, R. & Goupil, M.-J. 2001, A&A(in press) 
Samadi, R., Goupil, M.-J., Lebreton, Y., & Baglin, A. 

2000, in SOHO 10 / GONG 2000 Workshop (in press) 
Samadi, R., Houdek, G., Goupil, M.-J., & Lebreton, Y. 

2001, submitted to A&A 

Spiegel, E. 1962, J. Geophys. Res., 67, 3063 
Stein, F. & Nordlund, A. 1991, in Challenges to Theories 
of the Structure of Moderate-Mass Stars, Vol. 195 (ed 
D. Gough & J. Toomre (Springer- Verlag)) 
Stein, R. F. 1967, Solar Physics, 2, 385 



